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Quasiparticle states near the surface of a quasi-one-dimensional organic superconductor 
(TMTSF)2PF6 are studied based on an extended Hubbard model on a quasi-one dimensional lattice 
at quarter-filling. Three types of pairing symmetries, (i) p-wave, (ii) d-wave, or (iii) /-wave are as- 
sumed. The resulting surface density of states has characteristic features for each pairing symmetry: 
i) a zero-energy peak (ZEP) in a U-shaped structure, ii) a V-shaped structure without ZEP, and 
iii) a ZEP in a V-shaped structure. From these results, we propose that the tunneling spectroscopy 
serves as a strong method to identify the the pairing symmetry in (TMTSF)2PF6. 
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In recent years, pairing symmetry in various uncon- 
ventional superconductors, such as the high- Tc cuprates, 
heavy fermion systems, Sr2Ru04, organic supercon- 
ductors, and so on, has been exte*isively studied 

I both experimentally and theoretically.ETH In particu- 
lar, quasi-one-dimensional (QID) organic superconduc- 

; tors (TMTSF)2X (X=PF6, CIO4, etc.) have recently 
attracted much attention as a possible spin-triplet su- 
perconductor. Experimentally, the observation of a large 
critical magnetic field Hc2 exceeding Pauli paramagnetLc. 
limitia, as well as an unchanged Knight shift across T^B 

[ strongly suggest spin-triplet pairing. As for the orbital 
part of the pair wave function, the presence of nodes in 

^ the pair potential on the Fermi surface has been suga. 

'. gested from NMR measurements for (TMTSF)2C104lI 
and (TMTSF)2PF6,eI which exhibit the absence of Hebel- 
Slichter peak as well as a power-law decay of T^^ below 
Tc- On the other hand, a thermal conductivity measure- 

I ment has suggested the absence of nodes on the Fermi 

• surface in (TMTSF)2C104.B 

[ Theoretically, several previous-studies have proposed 
a triplet p-wave pairing state,Bilil for which the nodes 

' of the pair potential do not intersect the Fermi surface. 
On the other hand, a spin-singlet d-wave-like pairing me- 
diated by spia pfluctuations has been proposed by sev- 
eral authors This is because superconductivity lies 
right next to the 2kp spin density wave (SI)W) phase 
in the pressure-temperature phase diagram.113 Moreover, 
one of the present authors have recently proposedEEl 
that triplet /-wave-like pairing may dominate over d- 
and p-w&ve in (TMTSF)2PF6 due to a combination of 



QID Fermi surface, coexistence of 2fcF SDW and 2kp 
charge density|-wa]Ve (CDW) suggested from diffuse X- 
ray scattering ,113113 and an anisotropy in the spin fluctu- 
ations. 

Thus, the situation is not settled either experimentally 
or theoretically. The purpose of the present study is to 
propose an experimental method to determine which one 
of the pairing symmetries is realized in (TMTSF)2PF6. 

Now, for the high-Tc cuprates, which has a singlet d- 
wave pair potential, it has been clarified that the internal 
phase causes a drastic interference effect in the quasi- 
particle states near surfaces or interfaces, enabling us to 
detect the sign change in the pair potential. Namely, a 
zero-energy bound state (ZES) at a (110) surface of a 
d-wave superconductor reflects the sign-change of the ef- 
fective pair potential in the-jprocess of the reflection of 
quasiparticle at the surfacalj. The formation of ZES re- 
sults in a peak in the surface density of states (SDOS) 
at the Fermi energy (zero-energy) and manifests itself as 
a so-called zero-bias conductance pea|k4ZiBCP) observed 
in scanning tunneling spectroscopy,t3cj which is consid- 
ered as a strong evidence for a sign change in the pair 
potential. 

Recently, Sengupta et al. haice proposed that the pair- 
ing symmetry in (TMTSF)2XE£l can be determined from 
the presence/absence of the ZES on the surface. Al- 
though their study points out an important aspect, their 
argument is mainly restricted to the absence/presence of 
ZEP, from which the p-wave and /-wave pairings cannot 
be distinguished. In fact, as we shall see, one has to look 
into the overall line shape of the SDOS to distinguish p- 
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and /-wave pairings. Since the detailed line shape of the 
SDOS is significantly influenced by the actual shape of 
the Fermi surface, we have to consider a more realistic 
lattice structure, in which the quasi one-dimensionality 
(warping) of the Fermi surface is taken into account. 

In order to meet this requirement, we consider an 
extended Hubbard model on a QID lattice at quarter- 
filling, |-epctending the previous study on a 2D square 
lattice-Ea We concentrate on (TMTSF)2PF6 because 
there is no complexity (unit cell doubling) due to an- 
ion ordering like in (TMTSF)2C104.eI Three types of 
physically plausible pairings (i) triplet 'p-wave' (ii) sin- 
glet 'd-wave' and (iii) triplet '/-wave' are studied. The 
spatial dependence of the pair potentials is determined 
self-consistently, and the SDOS is calculated using the 
self-consistently determined pair potentials. We propose 
from the calculation results that the quasiparticle tun- 
neling spectroscopy should serve as a strong method to 
identify the pairing symmetry in (TMTSF)2PF6. 

The extended Hubbard model is given as 



H 



(iJ>b:Q 

(i,j)„,a,/3 i,Q 

where Cia [cj^] is the annihilation [creation] operator of 
an electron with spin a =],i at site i = (ia,ib)- Here 
ta[b] is the hopping integral, and (i, j)a[&] stands for sum- 
mation over nearest neighbor pairs in the a[6]-axis di- 
rection, respectively. V is the inter-electron potential 
between sites separated by m lattice spacings in the a- 
direction, and (i, j)m represents summation over pairs of 
sites separated by m lattice spacings. m depends on the 
choice of the pairing considered. We choose U/ta = 0.1 
in order to take into account the QID Fermi surface of 
(TMTSF)2PF6, which is open in the fcf,-direction. The 
chemical potential fi is determined so that the band is 
quarter-filled. 

By applying a mean-field approximation, A?^^ = 

■j(cictCj;3) is introduced, which represents the supercon- 
ducting pair potential for pairs formed by a-spin elec- 
tron on the ith site and /3-spin electron on the jth site. 
We assume that A^'' is proportional to Si^j^, where ib 
and jb are coordinates in the 6-direction. Thus the unit 
cell contains TVl sites in the a-direction and one site in 
the 6-direction. We consider three pairing symmetries 
shown in Fig. |l|. Namely, (i) p-wave: Sb = {Sb is the 
6-component of the total spin of a pair) triplet pairing 
between sites separated by 2 lattice spacings {rup = 2). 
This is a 'p-wave' pairing because the pair potential has 
a 2Ap sin 2kaa, fc-dependence in the bulk state, so that 
the pair potential changes its sign as H — along the Fermi 
surface (see Fig. ||). (ii) d-wave: singlet pairing between 



sites separated by 2 lattice spacings (m^ = 2). This 
is a 'c?-wave' pairing in the sense that the pair poten- 
tial changes its sign as H h— along the Fermi surface 

due to its 2 Ad cos 2kaa /c-dependence in the bulk state, 
(iii) /-wave: Sb = ±1 triplet pairing between sites sepa- 
rated by 4 lattice spacings {irif — A). This is an '/-wave' 
pairing in the sense that the pair potential changes its 

sign as H 1 1 — along the Fermi surface due to its 

2A / sin Akatt fc-dependence in the bulk state. These three 
pairings are physically plausible in the sense that they 
are consistent with the spin alignment of the 2/cf SDW 
phase of a quarter filled_^(fcF — j;^) system with an easy 
axis in the 6-direction.E3 In other words, the 2kp spin 
fluctuations can favor pairing in these channels. 

In order to represent A^J^ in a more convenient way, we 
introduce a new coordinate j along the a-direction. The 
original coordinate j is represented as j = (j,jfe) with j = 
1, • • • , A^L- In the 6-direction, we assume Nj, unit cells, 
and the electrons are Fourier transformed as Cja{kb) — 
E,1LlCjae-"=^^"^^ and C,p{-k,) = E,?=i Cj;3e^'=''^''^ 
where —ir/a < ki, < 7r/a, fc^ = j^n with n being an 
integer. 
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(i) p-wave 



(ii) d-wave 



(iii) f-wave 
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FIG. 1. (a) An illustration of Cooper pairing in real space 
and (b) the shape of Fermi surface for tb/ta = 0.1 at 
quarter-filling and the pair potential for (i) triplet 'p-wave' 
{nip — 2), (ii) singlet 'd-wave' (m^ = 2), and (iii) triplet 
'/-wave' {mf =4). In (a), the pairs are depicted by dashed 
lines, and in (b) -|- (— ) denotes the region where the sign of 
the pair potential is positive (negative). 



After the Fourier transformation, the mean-field 
Hamiltonian becomes 

hmf - J2 [^hC'b) clikb) a^i~kb) a^{~kb) 
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Here, for simplicity, the off diagonal part in Eq.(||) is 
assumed in the following forms. For the p-wave state. 
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at low temperatures, where pi^si'-^) denotes the SDOS at 
the i-th site from the surface in the superconducting state 
and /On(w) denotes the DOS in the normal state. In this 
paper, pn{uj) is obtained from the DOS at the Ni^ /2-th 
site far away from the surface. 



for the d-wave state, 
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and for the /-wave state, 

A^^ = A^^ = V A^.(5„±4, A^:^ = A^T = 0. 



(5) 
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We have taken the total number of sites as A^l — 500 
and Nb = 50. The value of the pair potential and 
the chemical potential in the bulk with V/ta — 4.0 are 
(i) Ap/ta = 0.280, fi/ta = 1.39, (ii) Aa/ta = 0.164, 
p./ta = 1.39, and (iii) Af/ta = 0.244, p./ta = 1.38, re- 
spectively. 

In the actual numerical calculation, the above 
Hamiltonian jt^vip is diagonalized by Bogoliubov 
transformationES given by Cia{kb) — X^iy ^j,i'7f (^fc)i ^^'^ 
Cjpi-kb) = Y.ulli^b)Ulj^^j,^, where v is the index 
which specifies the eigenstates. Then, the mean-field 
Hamiltonian described in Eq. (|^) is rewritten as Hmf = 
^kb.v Eiy{kb)^l{kb)^y{kb), where the operator 7^ (fcb) sat- 
isfies the fermion's anticommutation relation. The spa- 
tial dependence of the pair potential with ^-wave pairing 
symmetry is determined self-consistently as 



A' 



V 
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,{l-/[i?.(fc,)]}, (7) 



where j\Ei^{kb)\ denotes the Fermi distribution function. 
The procedure is iterated until the pair potential A'^- is 
obtained fully self-consistently. We calculate the SDOS 
using the pair potential determined self-consistently. In 
order to compare our theory with scanning tunneling mi- 
croscopy (STM) experiments, we assume that the STM 
tip is metallic with a flat density of states (DOS), and 
that the tunneling probability is finite only for the near- 
est site from the tip. This assumption has been verified 
through the study of tunneling conductance of unconven- 
tional superconductors. This is because the magnitude 
of the tunneling probability of an electron is sufficiently 
low in the actual STM experiments. The resulting tun- 
nelin&yonductance spectrum converges to the normalized 
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FIG. 2. The left panels are the spatial dependences of the 
pair potentials along the a-axis near the surface in the Z-wave 
state {I = p, d, /), and the right panels show the SDOS at the 
surface normal to the a- or 6-axis along with the bulk density 
of states. 



The obtained spatial dependences of the p-, d-, and 
/-wave pair potentials and the corresponding SDOS are 
plotted in Fig. ||. Since the spatial dependence of the pair 
potential is complex, we define the following quantities 
given by 



A^, 



Re[A^^^.±2]/Ap, 



A 



Re [A 
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0, 
0, 
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to visualize the spatial dependences clearly. The left pan- 
els of Fig. H is the obtained result for the spatial depen- 
dence of the pair potential near a surface normal to the 
a-axis, and the right panels show the SDOS at the sur- 
face normal to the a- or the &-axis along with the bulk 
DOS. 

First, let us look into the results for the triplet p- 
wave pairing state shown in the upper panels of Fig. ^ 
[see Fig. ||(a) and ||(b)]. Since triplet Cooper pair is 
formed between two electrons with 2 lattice spacings, 
A^^ — —A^^2,-a is satisfied. Both the magnitude of 
^ and A^ _^ is suppressed near the surface and ap- 
proaches 1 and —1 in the middle of superconductor, re- 
spectively. As shown in Fig. ||(b), the corresponding DOS 
has a U-shaped gap structure similar to that of the con- 
ventional s-wave pairing due to the fact that the nodes 
of the pair potential do not intersect the Fermi surface. 
The ZEP shows up in SDOS at the surface normal to 
a-axis due to the formation of ZES, since an injected 
and reflected quasiparticle feel different sign of the pair 
potentiaO. On the other hand, at the surface normal to 
the 6-axis, since an injected and reflected quasiparticle 
feel the same pair potential, the ZES is not formed and 
the resulting SDOS has no ZEP, resulting in a overall line 
shape similar to that of the bulk DOS. These results for 
the p-wave pairing are consistent with those in Ref. 

Next, we look into the corresponding quantities in sin- 
glet d-wave pairing case [see Fig. ||(c) and ||(d)]. Since 
the pair is formed between sites separated by 2 lattice 
spacings, A^^ = A^^2 -a satisfied. The obtained spa- 
tial dependence of the pair potential exhibits a atomic- 
scale spatial oscillation near the surface and converges to 
the bulk value toward the middle of the lattice. These 
features are similar to the previous results for—the ex- 
tended Hubbard model on a 2D square latticec^. The 
corresponding SDOS (and bulk DOS) has a V-shaped 
structure due to the existence of nodes of the pair poten- 
tial on the Fermi surface. However, since an injected and 
refiected quasiparticle feel the same pair potential both 
at the surfaces normal to a- and 6-axis, no ZEP appears 
at the SDOS [see Fig. ||(d)]. Due to the absence of the 
pair potential along the 6-axis, the ZEP never appears 
for arbitrary orientation of the surface which iSpStrikingly 
different from the case of the high- Tc cupratesE2l. 

Finally, we move on to the case of the triplet /-wave 
pairing. The /-wave pair [see Fig. ^(a) (iii)] is formed 
between sites separated by 4 lattice spacings, so the re- 



sulting pair potential satisfies A 



/ 



-A 



, . As seen 



from Fig. ||(e) , the obtained pair potential has a complex 
spatial dependence as compared to that of the p-wave 
pairing. Comparing Fig. ||(b) and Fig. ^(f ), it can be seen 
that the /-wave pairing belongs to the same class as that 
of the p-wave pairing as far as the absence/presence of 
the ZEP is concerned, as has been pointed out in Ref. |2^. 
However, since the /-wave pair potential has nodes on the 
Fermi surface, the resulting SDOS (and bulk DOS) has 
a V-shaped structure similar to that for the d-wave case 



in sharp contrast with the case of p-wave pairing. 

In total, as summarized in Table ^, the p-,d-,and f- 
wave pairings can be clearly distinguished from the com- 
bination of the overall line shape of the SDOS and the 
presence/absence of the a-axis ZEP. 

In summary, we have studied the quasiparticle SDOS 
of an organic superconductor (TMTSF)2PF6 based on 
an extended Hubbard model on a QID lattice at quar- 
ter filling. The non-local feature of the pair potential 
and the atomic-scale geometry of the surface are explic- 
itly taken into account in the present calculation. Three 
types of pairing symmetries, (i) p-wave, (ii) cZ-wave, and 
(iii) /-wave have been considered. The calculation re- 
sults suggest that we can clearly distinguish the present 
three pairing symmetries from tunneling spectroscopy. 
We believe our theoretical prediction can be verified ex- 
perimentally in the near future. It is an interesting future 
problem to investigate how our results will be modified 
for (TMTSF)2C104, in whicii a unit cell doubling due to 
anion ordering takes placeEH It is also an interesting fu- 
ture problem to study Josephson effect in QID organic 
superconductors since it has been clarified in the previ- 
ous studies that Josephson effect ace crucially infiuenced 
by unconventional pair potentials H 
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TABLE I. Surface density of states (SDOS) for p,d, and 
/-wave pairings. 



Symmetry 


SDOS 


p-wave 


U-shaped + a-axis ZEP 


d-wave 


V-shaped + No ZEP 


/-wave 


V-shaped + a-axis ZEP 
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